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Abstract 

We discuss quantum correlations in systems of indistinguishable particles in relation 
to entanglement in composite quantum systems consisting of well separated subsystems. 
Our studies are motivated by recent experiments and theoretical investigations on quan- 
tum dots and neutral atoms in microtraps as tools for quantum information process- 
ing. We present analogies between distinguishable particles, bosons and fermions in low- 
dimensional Hilbert spaces. We introduce the notion of Slater rank for pure states of pairs 
of fermions and bosons in analogy to the Schmidt rank for pairs of distinguishable parti- 
cles. This concept is generalized to mixed states and provides a correlation measure for 
indistinguishable particles. Then we generalize these notions to pure fermionic and bosonic 
states in higher-dimensional Hilbert spaces and also to the multi-particle case. We review 
the results on quantum correlations in mixed fermionic states and discuss the concept of 
fermionic Slater witnesses. Then the theory of quantum correlations in mixed bosonic 
states and of bosonic Slater witnesses is formulated. In both cases we provide methods 
of constructing optimal Slater witnesses that detect the degree of quantum correlations in 
mixed fermionic and bosonic states. 
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1 Introduction 



The understanding and characterization of quantum entanglement is one of the most funda- 
mental issues of modern quantum theory and a lot of work has been devoted to this 
topic in the recent years |3|-j6|. 

In the beginning of modern quantum theory, the notion of entanglement was first noted 
by Einstein, Podolsky, and Rosen Q, and by Schrodinger ||. While in those days quantum 
entanglement and its predicted physical consequences were (at least partially) considered as 
an unphysical property of the formalism (a "paradox"), the modern perspective on this issue is 
very different. Nowadays quantum entanglement is seen as an experimentally verified property 
of nature, that provides a resource for a vast variety of novel phenomena and concepts such as 
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quantum computation, quantum cryptography, or quantum teleportation. Accordingly there 
are several motivations to study the entanglement of quantum states: 



I. Interpretational and philosophical motivation: Entanglement plays an essential role in 
apparent "paradoxes" and counter-intuitive consequences of quantum mechanics |0-1O|. 



II. Fundamental physical motivation: The characterization of entanglement is one of the 
most fundamental open problems of quantum mechanics. It should answer the question 
what is the nature of quantum correlations in composite systems . 

III. Applied physical motivation: Entanglement plays an essential role in applications of 
quantum mechanics to quantum information processing, and in particular to quantum 
computing |H], quantum cryptography O] and quantum communication [p^] (i.e. 
teleportation [l^,[l6|] and super dense coding [jlj]). The resources needed to implement 
a particular protocol of quantum information processing are closely linked to the entan- 
glement properties of the states used in the protocol. In particular, entanglement lies at 
the heart of quantum computing [Q] . 

IV. Fundamental mathematical motivation: The entanglement problem is directly related 
to one of the most challenging open problems of linear algebra and functional analy- 
sis: Characterization and classification of positive maps on C* algebras (for 
mathematical literature see |21-p4[). 



While entanglement plays an essential role in quantum communication between parties 
separated by macroscopic distances, the characterization of quantum correlations at short 
distances is also an open problem, which has received much less attention so far. In this 
case the indistinguishable character of the particles involved (electrons, photons,...) has to be 
taken into account. In his classic book, Peres [|J discussed the entanglement in elementary 
states of indistinguishable particles. These are symmetrizations and antisymmetrizations of 
product states for bosons and fermions, respectively. It is easy to see that all such states of 
two-fermion systems, and as well all such states formed by two non-collinear single-particle 
states in two-boson systems, are necessarily entangled in the usual sense. However, in the case 
of particles far apart from each other, this type of entanglement is not of physical relevance: 
"No quantum prediction, referring to an atom located in our laboratory, is affected by the 
mere presence of similar atoms in remote parts of the universe" [|| . This kind of entanglement 
between indistinguishable particles being far apart from each other is not the subject of this 
paper. Our aim here is rather to classify and characterize the quantum correlations between 
indistinguishable particles at short distances. We discuss below why this problem is relevant 
for quantum information processing in various physical systems. Perhaps the first attempt 
to study such quantum correlations in macroscopic systems was done by A.J. Legett p5|. 



More recently he has formulated the concept of dis connectivity [26] of quantum states which 
is somewhat related to the concepts developed in this review. 

This paper is organized as follows: In section || we illustrate the consequences of indistin- 
guishability. In section |3] we describe analogies between quantum correlations in systems of 
two indistinguishable fermions, bosons, and two distinguishable parties where we concentrate 
on the lowest-dimensional Hilbert spaces that allow for non-trivial correlation effects. We 
derive the fermion and bosons analogues of recent results by Wootters 27], by Kraus and 
Cirac [E3] and by Khaneja and Glaser [S3]. Our results shed new light on a question posed 
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recently by Vollbrecht and Werner: Why two qubits are special [31| (see also [|32|). In section ^| 
we report further results on quantum correlations in pure states of indistinguishable fermions 
in higher-dimensional cases. Results on mixed fermionic states are summarized in section |5], 
and in section |6| we report further results on identical bosons. We conclude in section 0. 



2 Quantum correlations and entanglement 

2.1 Physical systems: Quantum dots and neutral atoms in microtraps 



Semiconductor quantum dots |33[] are a promising approach to the physical realization of 
quantum computers. In these devices charge carriers (e.g. electrons) are confined in all three 
spatial dimensions. Their electronic spectrum consists of discrete energy levels since the 
confinement is of the order of the Fermi wavelength. It is experimentally possible to control 
the number of electrons in a such a dot starting from zero (e.g. in a GaAs heterostructure |34}| ). 

When one wants to use quantum dots for quantum computation it is necessary to define 
how the qubit (i.e. the basic unit of information) should be physically realized. E.g. the 
orbital electronic degrees of freedom or the electron spin S can be chosen to form the qubit. 
An advantage of the latter approach is that the decoherence time is much longer for the spin 
than for the orbital degree of freedom (usually three orders of magnitude [33. 35]). 



The implementation of quantum algorithms needs single qubit and two qubit quantum 
gates |36l]. For the spin degree of freedom the former can be achieved by the application 



of a magnetic field exclusively to a single spin [37]. It is well-known [38] that arbitrary 



computations can be done if, apart from single qubit rotations, a mechanism by which two 
qubits can be entangled is available (the entangling gate VSWAP , together with single qubit 
rotations, can be used to produce the fundamental controlled-NOT gate). It was proposed 
in [37] to realize this mechanism by temporarily coupling two spins in two dots. The coupling, 
described by a Heisenberg Hamiltonian H{t) = J(t)S\S2, can be turned on and off by lowering 
and raising the tunnel barrier between neighboring quantum dots. 

Another interesting type of physical implementation possibilities are neutral atoms in 
magnetic |3J| or optical [4C] microtraps. Here each single neutral atoms is trapped in a 



harmonic potential and their collisional interaction can be controlled by temporarily decreasing 
the distance of the traps or by state-selective switching of the trapping potential. 



2.2 Consequences of indistinguishability 

We will use a schematic view of two electrons located in a double-well potential to illustrate 
the consequences of indistinguishability for entanglement. This description applies to the 
discussed examples of quantum information processing in quantum dots and in optical or 
magnetical microtraps (replacing electrons by atoms). For this illustration we will assume 
that the qubit is modeled by the spin degree of freedom, which we will denote by | f) and 
| {). Furthermore we have two spatial wavefunctions labeled and initially localized 
in the left and in the right potential well, respectively. Then the complete state-space is four 
dimensional: {|0)| |>, |x)l T), [0)1 I), lx)l !)}• 

We start with a situation where we have one electron in each well. Even if they are prepared 
completely independently, their pure quantum state has to be written in terms of Slater 
determinants in order to respect the indistinguishability. Operator matrix elements between 
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such single Slater determinants contain terms due to the antisymmetrization of coordinates 
("exchange contributions" in the language of Hartree-Fock theory). However, if the moduli 
of (r\<f>), (r\x) have only vanishingly small overlap, these exchange correlations will also tend 
to zero for any physically meaningful operator. This situation is generically realized if the 
supports of the single-particle wavefunctions are essentially centered around locations being 
sufficiently apart from each other, or the particles are separated by a sufficiently large energy 
barrier. In this case the antisymmetrization has no physical effect and for all practical purposes 
it can be neglected. 

Such observations clearly justify the treatment of indistinguishable particles separated 
by macroscopic distances as effectively distinguishable objects. So far, research in quantum 
information theory has concentrated on this case, where the exchange statistics of particles 
forming quantum registers could be neglected, or was not specified at all. 

Under these conditions we write an initial state \ipmit)AB = 1 0)1 ])a ® \ x) I V)B where A 
(Alice) and B (Bob) are (physical meaningful) labels for the particle in the left and the right 
dot, respectively. The situation is shown in figure [j]. 




Figure 1: The initial state \^imt)AB = \4>)\ 1)a (8> \x) \ 1)b i-e. the case of separated wells. \<j>) 
and \x) denote the spatial part of the wave- function localized in the left and in the right well, 
respectively. 

Now we want to analyze the situation when the two wells have been moved closer together 
or the energy barrier has been lowered. In such a situation the probability of finding, e.g., 
Alice's electron in the right well is non-vanishing. Then the fermionic statistics is clearly 
essential and the two-electron wave-function has to be antisymmetrized and reads \ip{t\)) = 
75 [ 10)1 T)i ® \x) \ ih - |x)| Di ® 10)1 D2 ]■ The indices A and B are changed to 1 and 2 
here to stress that the enumeration of the particles is completely arbitrary since these labels 
are not physical: because of the spatial overlap of the wavefunctions the individual particles 
labeled '1' or '2' are not accessible independently. The situation is shown in figure ^. 

Note that not only the labeling of the particles but also the notation suggesting a ten- 
sor product structure of the space of states is misleading because the actual state space is 



just a subspace of the complete tensor product [41]. As a consequence of this fact the anti- 
symmetrized state \ip(ti)) formally resembles an entangled state although it is clear that this 
entanglement is not accessible, and therefore cannot be used as a resource in the sense dis- 
cussed above for distinguishable particles. To emphasize this fundamental difference between 
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distinguishable and indistinguishable particles, we will use the term quantum correlations to 
characterize useful correlations in systems of indistinguishable particles as opposed to corre- 
lations arising purely from their statistics (thereby following p3fl). 

Quantum correlations in systems of indistinguishable fermions arise if more than one Slater 
determinant is involved, i.e. if there is no single-particle basis such that a given state of iV 
indistinguishable fermions can be represented as an elementary Slater determinant (i.e. a 
fully antisymmetric combination of N orthogonal single-particle states). These correlations 
are the analogue of quantum entanglement in separated systems and are essential for quantum 
information processing in non-separated systems. 

As an example suppose it is possible to control the coupling J(t) of the electrons such that 
at time t% 



m 2 )) = \ 



t)i ® ix)i i) 2 - ix)i 1)1 ® m t> 2 + m 1)1 ® ix)i t> 2 - ix)i t)i ® 1^)1 1) 2 ; 



which is illustrated in figure || 

In the given single-particle basis, {ipfo)} is written in terms of two elementary Slater 
determinants (and evidently there is no basis in which it can be written as a single one). This 
state contains some useful correlations beyond the required permutation symmetry as can be 
seen through localizing the particles again by switching off the interaction, i.e. raising the 
tunneling barrier or moving the wells apart (here we neglect the effects of non-adiabaticity, 



sec 



33fl for a more detailed study of these effects) . This corresponds to a partition of the basis 



between Alice and Bob, such that Alices Hilbert space is formed by {|</>)| f), \4>)\ j)} and Bobs 
by {|x)l t)> lx)l I)}- Then again the electrons can be viewed as effectively distinguishable, 
provided that none of the dots is occupied by two electrons. This does not happen here because 
the final final state is |Vwi)ab = 775 [ 10)1 X) A ® |x)| 1)b + \4>)\ Da ® |x)l T)b ], where new 
labels A and B are attributed to the particles, corresponding to the dot in which they are 
found after separation, i.e. the electron found in the left (right) dot is named A (B). \ip& n a,i)AB 
is shown in figure ^| 

The final state |^fi na l)AB> shown in figure [| is the Bell state \^ + ), i.e. a maximally en- 
tangled two qubit state (thus the operation performed in this example is the entangling gate 
VSWAP). In this sense it is reasonable to call {ipfe)) a maximally correlated state of two in- 
distinguishable fermions in a four-dimensional single-particle space and to view it as a resource 
for the production of entangled states of distinguishable particles. 



Motivated by these considerations in [43| we have developed a classification of states of 





Figure 2: Illustration of \ip(ti)}, i.e. after lowering the tunnel barrier. The electronic wave- 
functions are no longer completely localized in one of the wells. 
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Figure 3: Illustration of IVK^)) 




Figure 4: The final state |V>finai)AB- Raising the tunneling barrier localizes the wavefunctions 
again. 
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two fermions with 2M accessible single-particle states. This question was also addressed very 
recently by Li et al. [44] and by Paskauskas and You [45|. In these papers two-boson systems 
are examined and analogues to earlier results about two-fermion systems [ 33 , 43 1 are derived. 
However, Refs. [44] and differ in detail about which two-boson states should be considered 
as analogues of entangled states (in a bipartite system) in a certain limiting case. 

Zanardi [41,42] discusses another approach, ignoring the original tensor product structure 
through a partition of the physical space into subsystems. The entangled entities then are 
no longer particles but modes. This approach may be seen as complementary to the one 
followed here. For completeness we present the corresponding formalism in appendix |A| It is 
reasonable to consider both kinds of quantum correlations - which one is more useful depends 
on the particular situation, for instance on their usefulness for concrete applications, e.g. 
cryptography or teleportation. 



3 Analogies between bosons, fermions, and distinguishable par- 
ticles in low-dimensional Hilbert spaces 

3.1 Pure states: Schmidt rank and Slater rank 
3.1.1 Schmidt rank of distinguishable particles 

The "classic" examples for quantum entanglement were studied in systems composed of sep- 
arated (and therefore distinguishable) subsystems. The most investigated case involves two 
parties, say A(lice) and B(ob), having a finite-dimensional Hilbert space TLa and TLb, respec- 
tively. This results in a total space TL = Ha ®1~Lb- An important tool for the investigation of 
such bipartite systems is the bi-orthogonal Schmidt decomposition [fflj. It states that for any 
state vector \tp) G 7i there exist bases of TLa and TLb such that 

r r 

IV>> = ^2z i (\a i )®\bi)) 1 ^ > and J^z 2 = 1, (1) 

i=l i=l 

where the basis states fulfill that {ai\aj) = {bi\bj) = Sij. Thus, each vector in both bases 
for Ha and TLb occurs at most in only one product vector in the above expansion. The 
expression ([[]) is an expansion of the state into a basis of orthogonal product vectors 
with a minimum number r of nonzero terms. This number can take values between one and 
min{dim?^, dim 7^} and is called the Schmidt rank of \tp). \tp) is entangled if and only if 
r > 1. 



3.1.2 Slater rank of fermionic states 

Let us now turn to the case of two identical fermions sharing an n-dimensional single-particle 
space TL n . The total Hilbert space is A(TL n <g> TL n ) where A denotes the antisymmetrization 
operator. A general state vector can be written as 

n 

H = E <%///>> ( 2 ) 
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with fermionic creation operators /• acting on the vacuum |0). The antisymmetric coefficient 

-I (3) 



matrix Wij fulfills the normalization condition 



tr (w*w) = — — . 
Under a unitary transformation of the single-particle space 

J-f ] 
J jiJj 



w transforms as 



w 



UwU T 



(4) 



(5) 



where U T is the transpose of U . For any complex antisymmetric n x n matrix w there is a 
unitary transformation U such that w' = UwU T has nonzero entries only in 2 x 2 blocks along 
the diagonal [43, 4^]. That is, 



w 



diag [Zi,...,Z r ,Z ] with Z { 



-z, 



Zi 





(6) 



where z% > for i £ {1, . . . , r}, and Zq is the (n — 2r) x (n — 2r) null matrix. Each 2x2 block 
Zi corresponds to an elementary Slater determinant. Such elementary Slater determinants are 
the analogues of product states in systems consisting of distinguishable parties. Thus, when 
expressed in such a basis, the state \w) is a sum of elementary Slater determinants where 
each single-particle basis state occurs at most in one term. This property is analogous to the 
bi-orthogonality of the Schmidt decomposition discussed above. The matrix (||) represents an 
expansion of \w) into a basis of elementary Slater determinants with a minimum number r of 
non-vanishing terms. This number is analogous to the Schmidt rank for the distinguishable 
case. Therefore we shall call it the fermionic Slater rank of \w) |43| , and an expansion of the 
above form a Slater decomposition of \uu). 



3.1.3 Slater rank of bosonic states 

Similarly, a general state of a system of two indistinguishable bosons in an n-dimensional 
single-particle space reads 

n 

\ V ) = J2 Vl&}|0> (7) 

with bosonic creation operators b\ acting on the vacuum state. The symmetric coefficient 
matrix transforms under single-particle transformations just the same as in the fermionic 
case, 

v i — * UvU T . (8) 

For any complex symmetric matrix v there exists a unitary transformation U such that the 
resulting matrix v' = UvlI T is diagonal f44|-f46[| , i.e. 

v' = diag(^i, . . . ,z r ,0, . . . ,0) (9) 



9 



with Zi 7^ for i E {1, . . . , r}. In such a single particle basis the state \ v) is a linear combination 
of elementary two-boson Slater permanents representing doubly occupied states. Moreover 
Eq. (^) defines an expansion of the given state \v) into Slater permanents representing doubly 
occupied states with the smallest possible number r of nonzero terms. We shall call this 
number the bosonic Slater rank of \v). 

An expansion of the form @ for a two-boson system was also obtained very recently in 
Refs. [44,4^]. Moreover the fermionic analogue @ of the bi-orthogonal Schmidt decomposition 



of bipartite systems was also used earlier in studies of electron correlations in Rydberg atoms 



Regarding Slater determinant states in fermionic Hilbert spaces we also mention inter- 
esting earlier work by Rombouts and Heyde Q] who investigated the question under what 
circumstances a given many-fermion wavefunction can be cast as a Slater determinant built 
up from in general non- orthogonal single-particle states. Since a wavefunction of this kind can 
in general not be written as a single Slater determinant constructed from orthogonal single- 
particle states, i.e. has non-trivial quantum correlations beyond simple antisymmetrization 



effects, the criteria obtained in [48] do not address the issues here. 

As we saw, elementary Slater determinants in two-fermion systems, i.e. states with Slater 
rank one, are the natural analogues of product states in systems of distinguishable parties. 
One needs at least a fermionic state of slater rank two to form a quantum correlated state 
corresponding to a Schmidt rank two state of separated particles. In contrast, in the bosonic 
case one needs at least a state of Slater rank four to perform the same task. This can already 
be seen from the fact that for distinguishable particles entangled states need at least a 2 x 2 
dimensional Hilbert space. 

In the following we shall refer to the Schmidt rank or Slater rank of a given pure state also 
as its quantum correlation rank. 



3.2 Magic bases, concurrence, and dualisation 

One of the most important issues in Quantum Information Theory is the qualification and 
quantification of entanglement between several subsystems in a given state p. For the case 
of two distinguishable parties, a useful measure of entanglement of pure states is the von 
Neumann-entropy of reduced density matrices constructed from the density matrix p = \ip)(ip\ 

E iH)) = ~ tr A {pa log 2 pa) = -tr B (pb log 2 Pb) , (10) 

where the reduced density matrices are obtained by tracing out one of the subsystems, pa = 
tr bP and vice versa. With the help of the bi-orthogonal Schmidt decomposition of \tp) one 
shows that both reduced density matrices have the same spectrum and therefore the same 
entropy, as stated in Eq. (|To|) . In particular, the Schmidt rank of equals the algebraic rank 
of the reduced density matrices. A pure state is non-entangled if and only if its reduced density 
matrices are again pure states, and it is maximally entangled if its reduced density matrices 
are "maximally mixed", i.e. if they have only one non-zero eigenvalue with a multiplicity of 
minldimT^A, dimHs}- 
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3.2.1 Two distinguishable particles 



The lowest-dimensional system of two distinguishable particles having non-trivial entangle- 
ment properties consists of just two qubits, i.e. dimW^ = dim 7^ = 2. For this system, the 
entanglement measure ( |io|) takes a particularly simple form if the state \ip) is expressed in the 
so-called magic basis, 

Ixi) = \ (I T!> - I IT}) 
lx 2 ) = -j= (I TT) + I ID) 
|X3) = ^(ItD + IID) 

1x4) = ^ (ITT) -I ID) (ii) 

i.e. \ip) = ^2 i= i <Xi\xi)> where an obvious notation has been used. With respect to this basis 
it holds MM 



where h(x) = — xlog 2 (x) — (1 — x) log 2 (l — x) is the binary entropy function and the "concur- 
rence" C(\ip)) is defined by C(\i/j)) = \ Ylt=i a i\- Thus, a state is fully entangled if and only if 
all its coefficients with respect to the magic basis have the same phase. What is furthermore 
"magic" about this basis is the fact that its elements are (pseudo-)eigenstates of the time 



reversal operator D 27]: 

V\Xi) = \Xi) (13) 



with 



V = {{ia\) ® {ial))K . (14) 



Here a\ B are Pauli matrices in the basis used in the construction of the magic basis ([TI|) , and 
/C is the operator of complex conjugation which acts on a product state of basis vectors as 
K\fj,v) = \fiv), where fj,, v £ {f, J.}, and on a general vector as K^2^ v ip IM/ \fJ,i') = Yjp, v ^tvl^)- 
These relations are part of the definition of /C. T> is invariant under arbitrary SU(2) trans- 
formations performed independently on the two subsystems, due to the relation Ua y U T = a y 
for any U G SU(2). It was pointed out in Ref. [^TJ that this property is particular to the 2x2 
matrix a y , and does not have a true analogue in higher dimensions. However, in the following 
we will encounter similar invariance relations in higher-dimensional spaces where the manifold 
of transformations is restricted to a physically motivated subgroup of the unitary group. 

Using the time reversal operator T>, the concurrence can be written as C{\ip)) = 
with = T>\ip). Moreover, since the entanglement measure E{\ip)) is a monotonous function 
of C{\ip)) with both functions ranging from zero to one, one can equally well use C(\i/j)) as a 



measure of entanglement [27,51], as we shall do in the following. 
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The definition of the magic basis (and as well of the complex conjugation operator K) 
refers explicitly to certain bases in the two subsystems. However, using the above invariance 
property, it is straightforward to show that switching to different local bases has only trivial 
effects without any physical significance. In particular, the concurrence is invariant under 
such operations. This can be seen from writing \ip) in the computational basis as l^) = 

v =\ i ii,v\P"v) . Then C(\ip)) = 4| det(V')| 2 and in this context it has been named "tangle" 
by Wootters f28fl. 



3.2.2 Two fermions 

Let us now turn to the case of two fermions. The lowest-dimensional system allowing a Slater 
rank larger than one has a four-dimensional single-particle space resulting in a six-dimensional 
two-particle Hilbert space. This case was analyzed first in Ref. [33] where a fermionic analogue 



of the two-qubit concurrence was found. It can be constructed in the following way: For a 
given state \ w) = Ylt j=i w ijflfj\ty defined by its coefficient matrix Wij, define the dual matrix 
Wij by 



2 

k,l=l 



\ £ e ijkl w* kl , (15) 



with e ijfci being the usual totally antisymmetric unit tensor. Then the concurrence C(\w)) can 
be defined as Ip2| 



C(\w)) = \(w\w)\ 



e^WijWki 

i,j,k,l=l 



■ (WUW34 + + WUW23) I • (16) 



Obviously, C(\w)) ranges from zero to one. Importantly it vanishes if and only if the state \w)) 
has the fermionic Slater rank one, i.e. is an elementary Slater determinant. This statement 
was proved first in Ref. [j33[| ; an alternative proof can be given using the Slater decomposition 
of \w) and observing that 

detw = (l(w\w)) . (17) 



This relation is just a special case of a general expression for the determinant of an antisym- 
metric (2K) x (2K) matrix w, 



detW= (wKi. ^ £ n - t2K w hi2 ■■■w i2K _ li2K ^ , (18) 

il,...,l2A' = l 

which can be proved readily by the same means. Note the formal analogy to the definition of 
the tangle of two qubits [^] . 

Moreover, it is straightforward to show that the concurrence is invariant under arbitrary 
SU(4) transformations in the single-particle space. More generally, for any two states \w\), 
\u)2) being subject to a single-particle transformation U, \wi) 1— > \ui) = U\wi), i £ {1,2}, it 
holds 

(u2\ui) = detU (u)2 \wi ) . (19) 
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In summary, the concurrence C(\w)) constitutes, in analogy to the case of two qubits, a mea- 
sure of "entanglement" for any two-fermion state \w). It ranges from zero for "non-entangled" 
states (having Slater rank one) to one for "fully entangled" states which are collinear with 
their duals. In this case the Slater decomposition of \w) consists of two elementary Slater 
determinants having the same weight. 

From the occurrence of the complex conjugation in the definition of the dual state (^) it 
is obvious that the dualisation of a two-fermion state is an anti-linear operation. In fact, as 
we will illustrate below, the dualisation of a state \w) again corresponds to time reversal. As 
another physical interpretation, it can be identified with a particle-hole-transformation, 

u P hfX h = h , u P h\o) = f!flflfl\o) , (20) 

along with a complex conjugation. The operator of dualisation T>, \w) = T>\w), can be written 
as 

V = -U ph K (21) 

where KL is the anti-linear operator of complex conjugation acting on the single-particle states 
and the fermionic vacuum as 

!Cf}lC = f} , JCf i JC = f i , /C|0) = |0). (22) 

As a further analogy to the two-qubit case one also has a magic basis, i.e. a basis of 
(pseudo-)eigenstates of the dualisation operator (|2l|): 



Ixi) 


1 

71 


{flfl 


+ flfl) 


|0) 


|X2> 


i 


{flfl 


+ flfl) 


|0> 


|X3) 


1 

V2 


/1/4 


+ flfl) 


|0) 


1X4) 


i 

71 


{flfl 


— flfl J 


|0) 


|X5> 


i 

71 


{flfl 


— flfl J 


|0) 


Ixe) 


i 

71 


/1/4 


-flfl) 


10) 



Expressed in this basis, the concurrence of a state \w) = Yli=i a i\Xi) j us t reads C(\w)) = 
I Y^=i a i\- Therefore, \w) has a concurrence of modulus one if and only if all its coefficients 
with respect to the magic basis have the same phase. 

Like in the two-qubit case, the definition of the complex conjugation operator /C and the 
magic basis refers to a certain choice of basis in the single-particle space. However, due to the 
invariance properties described above, this does not have physically significant consequences. 

A four-dimensional single-particle space can formally be viewed as the Hilbert space of a 
spin-|-object. Having two indistinguishable fermions in this system means coupling two spin- 
| states to total spin states that are antisymmetric under particle exchange. This leads to 
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the two multiplets with even total spin, i.e. 5 = 2 and S = 0. When the single-particle states 
labeled so far by i S {1, 2, 3, 4} are interpreted as states of a spin-| particle with S z = |— (i— 1) 
(h = 1) these multiplet states read explicitly: 



2,2) 
2,1) 

2,0) 



A7JI0) 
flfl\o) 




(flfl + flfl)\o) 



2,-1) 
2,-2) 



flfl\o) 
flfl\o) 



0,0) 



(ftfl - M) io) 



(24) 



The phase of the last singlet state has been adjusted such that the dualisation operator reads 
in this ordered basis 



which is nothing but the time reversal operator acting on the two multiplets S = 2 and 
S = 0. Thus, the operation of dualisation can also be interpreted as a time reversal operation 
performed on appropriate spin objects. However, this interpretation is not necessarily the 
physically most natural one: the notion of dualisation and entanglement-like quantum corre- 
lations between indistinguishable fermions was first investigated for the case of two electrons 
(carrying a spin of ^) in a system of two quantum dots [|33|]. There the interpretation of T> 
as the anti-linear implementation of a particle-hole-transformation, rather than time reversal 
of formal spin objects, seems clearly more appropriate. Therefore we shall retain the general 
term dualisation for such an operation in the two kinds of systems investigated so far, and 
also in the case of a bosonic system to be explored below. 

3.2.3 Two bosons 

Having established all these analogies between two qubits and a system of two fermions sharing 
a four-dimensional single particle space, it is natural to ask whether similar observations can 
be made for a system of two bosons. The bosonic system showing properties analogous to 
those discussed before consists of two indistinguishable bosons in a two-dimensional single- 
particle space. This is the smallest-dimensional system admitting states with a bosonic Slater 
rank greater than one. It can be viewed as the symmetrized version of the two-qubit system, 
and its two-boson space represents therefore the Hilbert space of a spin-l-object. 

A general state vector of this system reads \v) = Ylt j=i v ijb\bj\0) with a coefficient matrix 



V = 



/ 1 o \ 

0-100 

1 
0-10 

1 

\ 1 / 



(25) 




(26) 
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being subject to the normalization condition 2|^4| 2 +4|i3| 2 + 2|C| 2 = 1. The appropriate 
dualisation operator is 

V = 1ZK (27) 

where K, is the complex conjugation operator with analogous properties as above, and the 
operator 1Z acts in the single-particle space as 

b\ -> Kb{rt = £ -ia\ a b\ , (28) 
6=1 

i.e. 1Z exchanges the state labels inferring a sign. When expressed in the ordered basis 
( (1/^2)44 1°)> 44 1°), (l/v/2)44l°) ) the dualisation operator reads 



V = -1 NC (29) 




which is, as natural, just the time reversal operator for a spin-l-object. Now the concurrence 
C(\v)) can be defined as the modulus of the scalar product of \v) with its dual \v) =V\v) and 
expressed as 

C(\v)) = \{v\v)\ = \4(AC-B 2 ) \ = |4detv| . (30) 

The last equation makes immediately obvious that C(\v)) is invariant under SU(2) transfor- 
mations in the single-particle space and is zero if and only if \v) has bosonic Slater rank one. 
Therefore this quantity has properties that are completely analogous to the two other cases 
before. A magic basis of the two-boson space is given by 

ixi) = ^ (44 + 44) |o> 
ix 2 ) = ~(44-44)|o) 

1X3) = i(b\b\) |0). (31) 



Again, the concurrence is maximal, i.e. C(\v)) = 1, if and only if the coefficients of \v ) in the 
magic basis have all the same phase. 

3.3 Mixed states and unified correlation measure: Wootters' formula and 
its analogues 

We now turn to the case of mixed states characterized by a density matrix p. For a general 
bipartite system (not necessarily consisting of just two qubits) the entanglement of formation 
of a state p is defined by the expression [|2, 54 1: 



£ F (p) = M {pi>im \ J2PiE(\i>i)) \ (32) 
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where E is the entanglement measure (|10|) for pure states, and the infimum is taken over all 
decompositions p = ')2iPi\'4 ) i){il ; i\ of the density matrix p in terms of normalized but in general 
not orthogonal states \ipi) and positive coefficients pi with YliPi = 1- 

In the case of two qubits an equivalent entanglement measure is, following Wootters [27|, 
given by 

C(p) = inf {p . W} |^p i C(|^})|. (33) 

where the concurrence C enters directly instead of via the entropy expression (^). This 
possibility relies on the fact that E(\ip)) for a pure state \ip) is a monotonous function of the 
modulus of its concurrence, and the infimum in ( j32[) is actually realized by a decomposition of 
p where each state has the same concurrence, such that the summation becomes trivial p7| . 

Furthermore C can be given by a closed expression avoiding the minimization entering 
(|33|), namely 

C(p) = max{0, Ai - A 2 - A 3 - A 4 } , (34) 

where Ai,A2,A3,A4 are, in descending order of magnitude, the square roots of the singular 
values of the matrix pp with p = T>pT>~ 1 . The singular values of the (in general non-hermitian) 
matrix pp can readily be shown to be always real and non-negative. 

The formula (p4h was first conjectured, and proved for a subclass of density matrices, by 



Hill and Wootters |50|| , and proved for the general case shortly afterwards by Wootters [27]. 

The expression (^3|) can obviously be read as a unified correlation measure for all of 
the three cases discussed above, namely two qubits, two indistinguishable fermions in a four- 
dimensional single-particle space, and finally two indistinguishable bosons in a two-dimensional 
single-particle space. This unified measure ranges from zero to one and vanishes for pure states 
with correlation rank one. Moreover this quantity is invariant under "local unitary transforma- 
tions", i.e. independent SU(2) transformation in the two qubit spaces, or, for indistinguishable 
particles, unitary transformations of the single-particle Hilbert space. Moreover, in all three 
systems, the dualisation operation has the property that for any two pure states \ip), \x) and 
there duals \x), the identity (V'lx) = (xl V 7 ) holds. 



These three properties are sufficient to take Wootters' proof [27| over to the other two cases. 



Therefore, the correlation measure C(p) of a general mixed state of two fermions sharing a 
four-dimensional single-particle space can be expressed as 

Cf(p) = max{0, Ai - A 2 - A 3 - A 4 - A 5 - A 6 } , (35) 

and the analogous expression for the two-boson system reads 

C B (p) = max{0, Ai - A 2 - A 3 } . (36) 

Again, in full analogy to the two-qubit case, the A, are, in descending oder of magnitude, the 
square roots of the singular values of of pp with p = VpV^ 1 being the appropriately dualised 



density matrix. As in Ref. [27 one can show that the singular values of pp are always real 
and non-negative. 
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In Ref. we have given a somewhat alternative proof of the Eqs. (|34D, (p5|), (36) in 
Ref. [43 1 by discussing the matrix 

Cij = (i>i\^j) (37) 

which arises from the decompositions p = \i/ji)(ipi\ of a given density matrix in terms of (in 
general non-normalized) states 



3.4 Invariance group of the dualisation and general unitary transformations 

To simplify the notation we shall in the following denote analogous quantities referring to the 
three systems above, namely two qubits, two fermions, and two bosons, in the form [•,-,•]. As 
an example, the unified correlation measure reads 

d 

C(p) =max{0,A 1 -^Ai}, (38) 

i=2 

where d = [4, 6, 3] is the dimension of the full Hilbert space of each system. 

Let us denote the group of (special) unitary transformations of the appropriate total Hilbert 
space by G = [SU(4), SU(6), SU(3)], and the group of "local transformations" (acting on the 
qubit spaces, or on the appropriate single-particle space) by H = [SU(2)<g)SU(2), SU(4), SU(2)]. 
In the following we do not distinguish between H and its [4, 6, 3] -dimensional representation 
in the full Hilbert space. 

For a given state |^) expressed in the magic basis, |^) = Yli=i a i\Xi)> the action of the 
dualisation operator is 

d 

V\-%) = Y,<\Xi) ■ (39) 

i=i 

Therefore, within the magic basis, the dualisation operator T> acting on the coordinate vector 
with components an is just V = l^xd^- We now investigate the subgroup M of G consisting 
of all unitary transformations IA that leave scalar products of the form (\I , i|\E'2) invariant, i.e. 
($i|$ 2 ) = (#i|* 2 ) for h-> = U\^i), i 6 {1,2} with two arbitrary states \^i). Such 
transformations must leave the dualisation operator invariant, 

UVU ] = V (40) 

which means in a magic basis 

UU T = 1 . (41) 

It follows that U is, in a magic basis, real and orthogonal, U G 0(d), d = [4,6,3], and, since 
detU = 1, U € SO(d). In a general basis U is an element of the corresponding equivalent 
representation of SO(d). Thus, M is just the special orthogonal group SO(d). Moreover, since 
H is certainly contained in M, and SO(d) is in all three cases a representation of H, we have 
M = H. This result can of course also be achieved by a straightforward construction of M in 
terms of its generators. 
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The fact that the group of "local transformations" is in all three cases identical to the 
unitary invariance group of the dualisation operator is peculiar to these small system studied 
here. As an example, consider a system of 2n qubits. The generalization of the operator T> is 
obvious, and, since the number of qubits is even, one also finds magic bases. Then the argument 
leading to the condition ( |40| ) can be repeated in a straightforward way, with the conclusion 
that the invariance group of T> is essentially SO(2 2n ), containing 2 2n_1 (2 2ri — 1) independent 
parameters. In contrast, H consists just of copies of SU(2) and has 6ra parameters. Thus only 
for n = 1 we find M = H; for larger n H is a proper subgroup of M. Similar arguments can 
be made for the generalization of the other two-particle systems to larger Hilbert spaces. 

As we saw, H contains in all of the three cases above all unitary transformations which 
turn the states of a given magic basis into another magic basis. This result is the key ingre- 
dient to prove a recent result by Kraus and Cirac for the two-qubit case [^] (also derived 
independently Khaneja and Glaser |]3(|). Here we shall just state the general result: Any 
unitary transformation U £ G of the full Hilbert space can be written in any magic basis as 

U = V x U d V 2 (42) 

with Vi,V2 € H and Ud € G diagonal. Thus, U d contains d — 1 = [3,5,2] independent 
parameters in the phase factors on its diagonal. Since the operations Vi, V 2 do not change the 
concurrence of any state, Eq. ([42]) constitutes a decomposition of a unitary transformation into 
parts which preserve "entanglement" and a part that does not. Therefore this decomposition 
is very helpful in the study of generation and deletion of entanglement in the course of a 



unitary transformation [29| 



The proof of the statement made in Eq. ([4^) for all three cases is completely analogous 
to the one given in the appendix of Ref. [^9| for the two-qubit case. For all systems the 
dimensions of the groups involved on the r.h.s. of Eq. (^) add up to give the dimension of 
the l.h.s. For analogous systems with larger Hilbert spaces such a decomposition of unitary 
transformations is in general not possible, as one can see from similar dimensional arguments 
as given above. 



3.5 Summary 

In this section we have studied three different two-body systems: 

• Two distinguishable qubits with a four-dimensional total Hilbert space. 

• Two indistinguishable fermions with a four-dimensional single-particle Hilbert space. 
The total two-particle space is six-dimensional. 

• Two indistinguishable bosons in a two-dimensional single-particle space. The total two- 
particle space is three-dimensional. 

These three systems share a series of properties: 

• They represent the smallest systems (regarding the dimensions of the Hilbert spaces) 
which allow for nontrivial correlation effects, i.e. a correlation rank larger than one. 

• In all three cases there is a correlation measure, the concurrence, related to an appropri- 
ately defined time reversal operator (dualisation). In the fermionic case this operator can 
be also interpreted as the anti-linear implementation of a particle-hole-transformation. 



18 



In all three cases a magic basis exists that consists of states being invariant under the 
equalisation operator. Maximally correlated states (with a concurrence of one) have real 
coefficients (up to a common phase factor) when expressed in the magic basis. 



• In the fermionic and bosonic case there are perfect analogues of Wootters' formula [27 
for the correlation measure of mixed two-qubit states. 

• In all three systems a general unitary transformation U. of the full Hilbert space can be 
written (in the magic basis) as 

U = ViU d V 2 (43) 

where Ud is diagonal and Vi,V2 are appropriate single-particle transformations. The 
result for the two-qubit case was derived recently in Refs. |p9| , |30|| . 

• In all three cases the group of single-particle transformations and the group of unitary 
invariance transformations of the dualisation operator coincide. In higher-dimensional 
cases the invariance group of the dualisation operator is larger than the group of single- 
particle transformations. 

Recently Verstraete et al. derived an explicit expression for a global unitary transformation 
which maximizes the entanglement of a given mixed state of two qubits Given the 



analogies derived above, we conjecture that the results of ref. |55|] can be applied to the other 
two systems studied here. 

In fact all these analogies can be derived from the following observation: as explained a 
two-fermion system in a four-dimensional single-particle space can be realized by two electrons 
on two neighboring quantum dots, taking into account only the lowest orbital states and the 
spin degree of freedom. After a projection onto the subspace where each dot contains one and 
only one electron a two-qubit system remains. Finally we arrive at a system of two bosons 
in a two-dimensional single-particle space by symmetrization. In this sense the system of two 
fermions is the most general one. Therefore, the question "Why two qubits are speciaF , posed 
recently by Vollbrecht and Werner |[l], can be formulated more generally: What is special 
about two fermions in a four-dimensional single-particle space? 



4 Quantum correlations of pure fermionic and bosonic states 
in higher-dimensional Hilbert spaces 

4.1 Two-fermion states 

As shown above, correlations of pure states of two fermions with a four-dimensional single- 
particle Hilbert space are characterized by the concurrence C(\w)). In particular this quantity 
vanishes if and only if \w) has Slater rank one. As soon as larger Hilbert spaces are considered 
there is no straight forward extension of the concurrence. However, performing the Slater 
decomposition is always possible and thus correlations of a two-fermion system can be charac- 
terized by the Slater rank. Here we introduce criteria for pure states in a single-particle space 
of arbitrary (even) dimension d = 2K that allow to specify the Slater rank without explicitly 
performing the decomposition. 
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To start with notice that in the case of 2K = 4 the concurrence vanishes if and only if 
the general expression ( |i~8|) for the determinant of an antisymmetric 2K x 2-ftT-matrix is zero. 
Using this for general dimension 2K together with the block structure of w in the canonical 
form it is clear that \w) has Slater rank < K if and only if 

2K 

^ ] E 1 2K Wi ^ _ _ ■ W i2K-l'i2K = 0- (44) 

Let us first turn to the case of 2K = 6. Here it only remains to distinguish Slater rank 
one from Slater rank two/three. This can be done by 

Lemma 4.1. Let \w) = j=i w ijfl fj\®) ^ e a two-fermion state in a six- dimensional single- 
particle space. Then \w) has Slater rank one if and only if Va, (3, with 1 < a < (3 < 6, 

6 

w ljWkl e^ Hal3 = 0. (45) 

i,j,k,l=l 



The proof is given in the appendix B.l. 



For 2K > 6 this result can be generalized to classify the Slater rank of an arbitrary pure 
state of two fermions. The quantities to be calculated have an analogous structure: Slater 
rank two or higher is distinguished from Slater rank one by contraction of two w-matrices 
with the e-tensor. Slater rank one and two can be distinguished from Slater rank three and 
above by contracting three ui-matrices and so on until finally states with maximal Slater rank 
are identified by the full contraction of K matrices. This is collected in the following lemma 
(see the appendix for the proof): 

Lemma 4.2. A two-fermion state \w) = Y^i^=i w ijfi fj\®) e ^^Hik ®7~L2k) has Slater rank 
< N if and only if Voi,... ,&2(k-N)> 1 < a.\ < . . . < a 2 (K-N) — 2K, 



2K 

E 

•• ,«2iV = l 



w i2N _ ll2N £ ll - l2Nai - a ^ K -^ =0. (46) 



4.2 Three- and multi-fermion states 

Three identical fermions in a 2iC-dimensional single particle space H.2K live in a total Hilbert 
space A(Ti.2K ® Ti-2K <8> T~(-2k)- In the usual notation a general pure state is written 

2K 

h = w ijkflfUl\°)> ( 47 ) 

i,j,k=l 

where w is completely antisymmetric and fulfills the normalization condition 

2K 

^2 wtjkWijk = g • (48) 

i,j,k=l 

Because for fermions every single-particle state can only be occupied once, a particle hole 
transformation (c.f. eqn. (^0|)) allows to treat the cases d = 2K = 4 and d = 5 more easily: 
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three fermions acting in a four-dimensional single-particle space are mapped by U p h onto one 
fermion (and thus such a state is never correlated) while a five-dimensional three-fermion state 
is mapped onto two fermions in five dimensions. This five-dimensional state can then always 
be embedded into a four-dimensional single particle space because the rank of the coefficient 
matrix is always even. Thus for d = 4 a state of three fermions is always non-correlated, for 



d = 5 the concurrence from section |3.2| can be used. The lowest non-trivial systems of three 
fermions act in a single particle space of dimension 2K = 6. 

In the case of distinguishable particles it is well known that a straight forward general- 
ization of the bi-orthogonal Schmidt decomposition to systems of more than two parties is 
not possible |5q-58|. Not every three-quantum-bit state can be transformed via local unitary 



transformations into a sum of tri-orthogonal product states, i.e. 

= do|0>A ® |0) B 8> \0)c + d!\l) A ® |1)b ® |l)c (49) 

is not always possible. 

Counting of dimensions shows that also for three fermions a canonical form consisting of 
a sum of elementary Slater determinants cannot exist. In the case of 2K = 6 this would mean 
to find a unitary transformation U of the single particle space such that 

2K 

™'lmn= E v>ijkUiiU mj U nk = Zl e lmn4 ™ + z 2 e 123lmn , (50) 

i ,j,k=l 

i.e. \w') oc z\ /1/2/3IO) + Z2 flflfl |0), where by adjusting the phases of the basis states the z k 
can be chosen real and non- negative. In general such a form depends on |_2-fT/3j real numbers 
([x\ is the largest integer < x), while a general three-fermion state has 



2K 
3 



(2iQ 3 -3(2iQ 2 + 4if 



parameters (including normalization and an overall phase) . A unitary transformation IA of the 
2/iT-dimensional single particle space has (2K) 2 real parameters. For the case of 2K = 6 this 
means that a general state is described by 40 parameters while the desired Slater decomposed 
form only has two parameters. But then the 36 real numbers that can be chosen for the unitary 
transformation are not enough in general to solve the corresponding set of linear equations. 

It is interesting that questions on multi- fermion pure states can be answered by projecting 
them onto pure states of fewer particles. For this purpose let us define an operator R\ a \ 
with \a) = ^2i a ifj\fy, which projects the three-fermion state \w) on a two-fermion state \w) 
supported at most on the single-particle space orthogonal to \a) and which is defined through 
the following equation: 

2K 2K 2K 

\w) = R\ a) \w) = R la) £ "•„,/;./•:/,() :>, ^ £w ijk a k f}f}\0). (52) 

i,j,k=l i,i=l k=l 

The first thing to be noticed is that a decomposition of \ w) into elementary non-overlapping 
Slater determinants only exists if there is a set {|e a )} Q=1 2 k-> (e Ql |e a2 ) = S aiCt2 and TL2K = 
span{|e a )} such that R\ e a\\w) has Slater rank one or is zero and R\ e a\\w) is orthogonal to 
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R\ e p\\w) if a ^ (3. In the basis where \w) has the desired block form this set is given by the 
set of basis vectors. 

Following the definition of non-correlated states for two-fermion systems a three-fermion 
state is non-correlated if and only if it can written as 

H = flflfl\0) (53) 

in some basis. This is simply a Slater decomposable state with Slater rank one. These are 
also the states called "separable" in reference |44]. We first remark that here the particle- 
hole transformation does not help to identify or even to quantify quantum correlations as 
it did via the concurrence in the two-fermion case, because for a three-fermion state \w) we 
have |(i/)|w)| = \(w\(UJC)^\w}\ = 0. This is analogous to three distinguishable qubits where 

MM) ® K?) ® (**c)W) = (cf. eqn. ©) 0. 

Instead of a direct translation of the concurrence there is a necessary and sufficient criterion 
for three-fermion pure states having Slater rank one, by using the projection defined above 



(the proof is given in appendix B.l 



Lemma 4.3. A three-fermion state \w) = Y^j,k=l w ijkfj f] /||0) G -4(W2K <8> H2K ® is 
noi correlated (has Slater rank one) if and only i/V|a) S 'H.2K'- 

2K 2K 
i,j=l k=l 

has Slater rank one or is equal to the 2K -dimensional zero. 



This lemma can be combined together with |4.2| (which allowed to calculate the Slater rank 
of two-fermion states) to give 

Lemma 4.4. A three-fermion pure state \w) = Y^l^ k=i w ijkfi fj fl\fy has Slater rank one if 
and only i/V«i,... ,(*2(K-2)> 1 < ai < • • • < &2{K-2) — 2if, and V a G C 2K , 

2K 

W ijk <H W lmn a t £ 0kmn ax ...a 2(K _ 2) = Q _ ^ 

i,j,k,l,m,n=l 

As an example consider 
k)=^/i t 4/3 t |0) + y/ 3 t / 5 t / 6 t |0), w l]k = l -(xe^ +ye 12 ^ k )\0}, |x| 2 + \y\ 2 = 1. (56) 



By lemma 4.3, \w) does not have Slater rank one because for a = (0, 0, 1,0, 0, 0) T 

R la) \w)xxfifi\0)+ y flfl\0) (57) 



has Slater rank two. Also evaluation of lemma [4.4| shows that \w) does not have Slater rank 
one because it gives rise to the five conditions 

a\xy = 0103x1/ = 0203x^ = 0305x17 = 0306 xy = (58) 
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which are obviously violated by some choices of a as long as x and y are non-zero. Notice 
that also \w) cannot be written by a transformation of the single-particle space as a sum of 
orthogonal Slater determinants as 

\w') oc Zl (f'iif'iif'ilO) + z 2 (/OlC/Ok/OllO), i.e. w' oc * /f/J/JlO) + ^2 44410), 

(59) 

because i?|^|u>') 7^ for all \b) but i2|^|io) = for \b) = /^|0) (or, equivalently, rank(w') > 
rank(w)). All the comments and results for three- fermion states can be carried over to systems 
consisting of N identical fermions in a Hilbert space A{Ti2K ® ■ ■ ■ ® T~L 2 k)- Here the lowest 
non-trivial dimension of the single-particle space is 2N since otherwise a treatment through a 
particle hole transformation is easier. A general pure state reads 



H= Yl w h-iNfh---fi N \°), Yl w *h-iN w ii-iN = Jfi> ( 60 ) 
«i,...,jjv=l ii,...,ijv=l 

where w is antisymmetric in all indices. Again for N > 2 a decomposition into orthogonal 
elementary Slater determinants is not possible in general, but a state equivalent to a pure 
separable state of distinguishable particles can be said to have Slater rank one if, by a trans- 
formation of the single-particle space, it can be cast into the form 

\w')= 4-../i|0). (61) 

This is a Slater decomposable state with Slater rank one and the definition again coincides 
with the definition of "separable" in reference \i4\ The operator R\ a \ can be generalized to a 
projector from a system of iV to N — 1 fermions as 

2K 

R \a)\w) = \w) with w il ... ijv _ 1 = N ^2 w h-i N ai N (62) 

ijv=l 

and it is clear that the necessary criterion for the existence of a decomposition into orthogonal 
Slater determinants given above for three fermions also holds for more particles. Furthermore 
lemma 4.3 can be translated as follows: 

Lemma 4.5. A pure N -fermion state has Slater rank one if and only if 

2K 

y\a)eH 2 K: R\ a) \w)=N Yl ^i..-<^4---/LJ°> ( 63 ) 

ii,... ,ijv=l 

has Slater rank one or is zero. 



The proof can be found in appendix BA. The lemma is maybe easier to handle if it is 
given in the following form: 

Lemma 4.6. A pure N -fermion state has Slater rank one if and only if V a, ; with 1 < a\ < 
■ ■ ■ < &2(K-2) < 2K, and for alia 1 ,... , a N ~ 2 G C 2K , 

2K 

£ w n ... lN ai . . . <"_ 2 2 a;,.., A a) x . . . af N Z 2 2 £ ^™-^ = . (64) 

«i,...ijv=l 
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Here lemma 4.5 is simply used several times, thereby reducing the number of fermions 
in the system to N — 1, N — 2 and so on until lemma [O^ can be applied to the resulting 
two-fermion state. 

As an example illustrating this lemma and showing that the straightforward extension of 
the two-fermion concurrence using the particle-hole transformation fails, let 

2K 

H= ^mf}f}flfl\0), w i3kl = x^™™ + ye^ kl ™ + ze 12 ^™ . (65) 

i,j,k,l=l 

To see that this is a correlated state take a = (0, 1, 0, 0, 0, 0, 0, 0) T to get R\ a )\w) = x f\f\f\\fy + 

y/|/|/g|0), which after some permutation of the basis is proportional to the three-fermion 
correlated state from equation (]56|). Thus also \ w) does not have Slater rank one. Nevertheless 

2K 

\w) = U ph K\w)= J2 c ijklmn ° P ™*ijMflflflfl\0) (66) 

i,j,k,l=l 
m,n,o,p=l 

and therefore (w\w) = although \w) contains correlations. 



4.3 Two-boson states 

Next we will apply the concept of assessing correlations by identifying the Slater number to 
bosonic pure states. Given a two-boson state \v) = Yl v ij°\°\ with Vij = vji a distinction 
between maximal and non-maximal Slater rank is provided by the vanishing of detv in the 
latter case, i.e. by considering the full contraction of K u-matrices with two i^-dimensional 
e-tensors. For a finer distinction in the case of a three-dimensional single-particle space we 
have 

^ o J. J. 

Lemma 4.1. Let \v) = 22i j=\ v ijK^j\^) ^ e a two-boson state in a three-dimensional single 
particle space. Then \v) has Slater rank one if and only if VI < a < 3 : 

3 

^ v ij v kl e lka ^ la = 0. (67) 

i,j,k,l=l 

This can be extended to higher-dimensional single-particle spaces essentially as for fermions. 
The result is the following lemma which is the bosonic equivalent to lemma |4.2| (the proof is 
sketched in appendix |B.2j ): 

Lemma 4.2. Let \v) be a two-boson state in A(TCk ®T~ik)- \v) has Slater rank < N if and 
only if Vai, . . . , uk-n, with 1 < a± < olk-n < K, 

K 

E„. . . . c .»it3."iajV-iai«JC-JV c .«2i4-..i2JV«iaif-iV _ n 

u %\%% U IZH ■ ■ ■ u l2N-l'l2N C C — U. V uo 7 

Hi- >*2JV = 1 
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4.4 Three- and multi-boson states 

For three bosons in a two-dimensional single-particle space a general state \v) has 8 real pa- 
rameters while a unitary transformation of the single-particle space has 4. This shows that, 
similar to the fermionic case, a straightforward generalization of the two-boson Slater decom- 
position is not possible, i.e. \v) cannot always be written as \v) = {z\b\b\b\ + Z2 62^2^2) 1 0) - 
The same holds for iV-boson states with N > 3. 

A special case of iV-boson states that indeed are Slater decomposable are non-correlated 
states that we are going to define now by generalizing the definition for two-boson systems: 

Definition 4.1. A pure N-boson state \v) is non- correlated (is Slater decomposable and has 
Slater rank one ) if in some basis it can be written 



L b\ . . . b\\0), i.e. v h ... iN = -L f[ Sijl . (69) 

It is worth noticing that this definition differs slightly from the one in reference where 
\v) is also called non-correlated if it can be written \v) = h= b\ . . . b N \0) in some basis where 
the bci are orthogonal. 

If we define an operator R\ a \ for \a) = J2i a «^l|0) that projects from a system of iV to a 
system of iV — 1 bosons (c.f. equations (|52|) and fl62])), 

K K 
\v)=R\a)\v)= Yl v n-iM^b J b\ i ...b\ N \0)=N %-^«iiv<---C-J ) ( ? °) 

i*JV)J=l ii...ijv=l 

(here the result is not necessarily only supported on the single-particle space orthogonal to 
I a)), then an identification of non-correlated states is possible. 

Lemma 4.3. A pure N-boson state \v) has Slater rank one if and only if the N — 1-boson 
state 

2K 

R\a)\v)=N ]T ^ 1 ... iw a ijv 6f l ...6t jv _ i | ) (71) 

il,... ,i]v=l 



has Slater rank one or is zero V|a) E H-2K- 



It is worth to notice that this lemma is also valid if, instead of the definition 4T , the 
definition from reference [44 is used. However, concentrating on definition |4.1| , it is possible 
to iterate lemma 4.3 until a two-boson state is left and lemma 4.2 can be applied. Thus we 
have 

Lemma 4.4. An N-boson pure state has Slater rank one if and only if Vctj with 1 < ct\ < 
■ ■■ < u K -2 < K, and for all a 1 ,... , a N ~ 2 £ C 2K , 

K 

V v h . . . a\ N ' 2) v h -a^ ...a"- 2 e iN - dN ~ iai - aK ~ 2 e iNjNai - aK -2 = 0. (72) 

«i,...ijv=l 
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5 Quantum correlations in mixed fermionic states 



In this section we repeat for completeness of the present paper the main results of Ref. [|43| 
concerning mixed states of two identical fermions. Part of these results, concerning the two- 
fermion system with a four-dimensional single particle space has been already discussed in 
Section 2. In this case the question of the Slater rank of a given mixed state can be answered 
exactly. Here we concentrate thus on the general case, for which no exact answers are known 
in general. 



5.1 Slater number of mixed fermionic states 

Let us now generalize the concepts introduced in Section 2 to the case of fermionic mixed states 
with a single-particle space of arbitrary dimension. To this end, we define formally the Slater 
number of a mixed state, in analogy to the Schmidt number for the case of distinguishable 



parties p9, 60] 



Definition 5.1. Consider a density matrix p of a two-fermion system, and all its possible 
convex decompositions in terms of pure states, i.e. p = YliPiWl* {^PlT > where denotes 
the Slater rank of |V ; )[ l ; the Slater number of p, k, is defined as k = min{r max }, where 
''max is the maximum Slater rank within a decomposition, and the minimum is taken over all 
decompositions. 

In other words, there exists a decomposition of p that uses pure states with Slater rank not 
exceeding k, and there is no construction that only uses pure states with Slater rank smaller 
than k. 

Many of the results concerning Schmidt numbers can be transferred directly to the Slater 
number. For instance, let us denote the whole space of density matrices in A{C 2K (g> C 2K ) by 
Six, and the set of density matrices that have Slater number k or less, by Slk- Slk is a convex 
compact subset of Six- A state from Slk will be called a state of (Slater) class k. Sets of 
increasing Slater number are embedded into each other, i.e. Six C SI2 C ...Slk... C Six- In 
particular, Sl± is the set of states that can be written as a convex combination of elementary 
Slater determinants; SI2 is the set of states of Slater number 2, i.e. those that require at least 
one pure state of Slater rank 2 for their formation, etc. 

The determination of the Slater number of a given state is in general a very difficult task. 
Similarly, however, as in the case of separability of mixed states of two qubits (i.e. states in 
C 2 <g> C 2 ), and one qubit-one qutrit (i.e. states in C 2 (g) C 3 ) the situation is particularly 
simple in the case of small K. For K = 1 there exists only one state (a singlet). For K = 2 we 
have presented in section 2 a necessary and sufficient condition for a given mixed state to have 
a Slater number of one. We have also defined the Slater correlation measure in analogy to 



Wootters' result |27|| , relating separability to a quantity called concurrence, and to eigenvalues 
of a certain matrix. One should note, however, that in the considered case of fermionic states 
there exists no simple analogy of the partial transposition, which is essential for the theory 



of entangled states of distinguishable particles. In fact, the Peres-Horodecki criterion [61,^2] 
in 2 x 2 and 2x3 spaces says that a state is separable if and only if its partial transpose is 
positive semidefinite. The results of section 2 can be formulated as the following Theorem: 

Theorem 5.1. Let the mixed state acting in A(C 4 (8> C 4 ) have a spectral decomposition 
p = Yui=i l^iK^il' where r is the rank of p, and the eigenvectors belonging to nonzero 
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eigenvalues Aj are normalized as = Aj5y. Let = ^ 6 w„f ) fafl |0) in some basis, 

and define the complex symmetric r x r matrix C by 

C y = (73) 

abed 

which can be represented using a unitary matrix as C = UCdU T , with = diag[ci, C2, ... , c r ] 
diagonal and |ci| > \c2\ > • • • > |c r |. The state p has Slater number 1 if and only if 

r 

|ci|<X)N- ( 74 ) 

The above Theorem is thus an analogue of the Peres-Horodecki-Wootters result for two- 
fermion systems having a single particle space of dimension 2K < 4. The situation is much 
more complicated, when we go to K > 2; this is similar to the case of the separability problem 
in C M <S> C N with M N > 6. These issues are investigated in the next subsection. 

At this point it is worth stressing that our formalism can also be used for odd dimensions, 
i.e. formally for half-integer K. Here the situation is again particularly simple for K = 1.5 
(2K = 3). All states have obviously the Slater number one and can thus be written as 
a convex sum of projectors onto states equivalent to the singlet state |0)|1) — |1)|0), which 
implies that the entanglement of formation of any antisymmetric mixed state (treated as a 
state of two distinguishable particles) is 1 ebit. This has been observed recently by Vidal, Diir 
and Cirac [p3[| , who showed that also the entanglement cost of every such antisymmetric state 
is 1 ebit which implies additivity of the entanglement of formation. 



5.2 Fermionic Slater witnesses 

We now investigate fermionic systems with single-particle Hilbert spaces of dimension 2K > 4. 
In this case, a full and explicit characterization of pure and mixed state quantum correlations, 
such as given above for the two-fermion system with K = 2, is apparently not possible. 
Therefore one has to formulate other methods to investigate the Slater number of a given 
state. We can, however, as in Ref. |13| follow here the lines of the papers that we have written 
on entanglement witnesses |34 , 65|, and Schmidt number witnesses [ 30 . In the following we 



skip the proofs, which are discussed in detail in Ref. [43|. 

In order to determine the Slater number of a fermionic density matrix p we note that due 
to the fact that the sets Slk are convex and compact, any density matrix of class k can be 
decomposed into a convex combination of a density matrix of class k — 1, and a remainder 



Proposition 5.1. Any state of class k, pk, can be written as a convex combination of a 
density matrix of class k — 1 and a so-called k— edge state 8: 

Pk = (1 -p)Pk-i +P$, 1>P>0, (75) 
where the edge state 5 has Slater number > k. 



The decomposition (|75|) is obtained by subtracting projectors onto pure states of Slater 
rank smaller than k, P = \il) <k )(ijj <k \ such that p^ — XP > 0. Here \t/j <k ) stands for pure 
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states of Slater rank r < k. Denoting by K(p), R(p), and r(p) the kernel, range, and rank of 
p, respectively, we observe that p' oc p — A|V')(V'I is non negative if and only if G R{p) and 
A < (^|p _1 |V>) _1 (see |6(|). The idea behind this decomposition is that the edge state 5 which 
has generically lower rank contains all the information concerning the Slater number k of the 
density matrix p^. 

As in the case of the Schmidt number, there is an optimal decomposition of the form ([75|) 
with p minimal. Alternatively, restricting ourselves to decompositions pk = SiPi!^)!'^!? 



with all rj < k, we can always find a decomposition of the form (75) with 5 £ Slk- We define 
below more precisely what an edge state is. 

Definition 5.2. A fermionic A:-edge state 5 is a state such that 5— e\ip <k )(ip <k \ is not positive, 
for any e > and \tp <k ). 

The edge states are thus characterized by the following 

Criterion 5.1. A mixed state 5 is a k-edge state if and only if there exists no \ip <k ) such 
that \ip <k ) e R{5). 

Now we are in the position of defining a fermionic A;-Slater witness (/c-SlW, k > 2): 

Definition 5.3. A hermitian operator W is a fermionic Slater witness (S1W) of class k if and 
only if Ti(Wa) > for all a 6 Slk-i, and there exists at least one fermionic p £ Slk such that 
Tr(Wp) < 0. 



It is straightforward to see that every S1W that detects p given by (75) also detects the 
edge state 5, since if TviW p) < then necessarily Tr(W5) < 0, too. Thus, the knowledge of 
all SlW's of fc-edge states fully characterizes all p S S7/%. Below, we show how to construct 
for any edge state a S1W which detects it. Most of the technical proofs used to construct and 



optimize fermionic Slater witnesses are practically equivalent to those presented in Ref. [34 
for entanglement witnesses. 

All the operators we consider in this subsection act in A{C 2K ® C 2K ). Let 5 be a k-edge 
state, C an arbitrary positive operator such that Tr(<5C) > 0, and P a positive operator whose 
range fulfills R{P) = K(5). We define e = inf|^< fe ^ (ip <k \P\ip <k ) and c = sup C\ip). Note 
that c > by construction and e > 0, because R(P) = K{5) and therefore, since R{5) does 
not contain any \if) <k ) by the definition of edge state, K(P) cannot contain any \4> <k ) either. 
This implies: 

Lemma 5.1. Given a fermionic fc-edge state 5, then 

W = P--C (76) 
c 

is a fermionic A>S1W which detects 5. 

The simplest choice of P and C consists in taking projections onto K(5) and the identity 
operator on the antisymmetric space l a , respectively. As we will see below, this choice provides 
us with a canonical form of a fermionic /c-SlW. 

Proposition 5.2. Any fermionic Slater witness can be written in the canonical form: 

W = W-el a , (77) 
such that RiW) = K(5), where 5 is a k-edge state and < e < mi^<k^(tp\W\ , ip) . 
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It is useful to consider also 



Definition 5.4. A fc-Slater witness W is tangent to Slk-i at p if 3 a state p £ Slk-i such 
that Tr(W» = 0. 

Similar methods to the ones used in ref. 
observe that 

Observation 5.1. The fermionic state p is of Slater class — 1 if and only if for all fermionic 
A;-SW's tangent to Sl k -i, Tr(W p) > 0. 



43 1 allow to prove the Proposition 5.2 and to 



5.2.1 Optimal fermionic Slater witnesses 



We will now discuss the optimization of a fermionic Slater witness. As proposed in [64| and [30 
an entanglement witness (Schmidt witness) W is optimal if there exists no other witness that 
detects more states than it. The same definition can be applied to fermionic Slater witnesses. 
We say that a fermionic k— Slater witness Wi is finer than a fermionic A:— Slater witness W\, 
if Wi detects more fermionic states than W\. Analogously, we define a k— Slater witness W 
to be optimal when there exists no finer witness than itself. Let us define the set of pure 
states of Slater rank k — 1 for which the expectation value of the fc-Slater witness W vanishes: 

T w = {|V>) such that (V>|W|V>) = 0} , (78) 

i.e. the set of pure tangent states of Slater rank < k. W is an optimal fc-SlW if and only 
if W — eP is not a &-S1W, for any positive operator P. If the set T\y spans the whole 
Hilbert space A(C 2K <g> C 2K ), then W is an optimal fc-SlW. If T w does not span A(C 2K <g> 
C 2K ), then we can optimize the witness by subtracting from it a positive operator P, such 
that PT\y = 0. For example, for Slater witnesses of class 2 this is possible, provided that 
inf| e)eC 2x [Pe 1/2 W e Pe 1/2 Uin > 0. Here for any X acting on A(C 2K ® C 2K ) we define 

l(e,.|A|e,.) - (e,.|X|.,e) 

- (., e\X\e, .) + (., e|A|., e)] , (79) 

as an operator acting in C 2K , and [A] m j n denotes its minimal eigenvalue (see |54|). An example 
of an optimal witness of Slater number k in A(C 2K (8) C 2K ) is given by 

W=la- j^P , (80) 

where V is a projector onto a "maximally correlated state", |^) = -7= YaLi HiHi+l^)- ^ ne 
reader can easily check that the above witness operator has mean value zero in the states 
/Ij/lj+ilO) f° r * = 1] 2, but also for all states of the form pjg^lO) where 

+f\e ilp21 + f\e l ^\ (81) 

g\ = _/t e -^12 +/ t e -^U 

-/jje-^+Zje-*** 1 , (82) 

for arbitrary ipij, i,j = 1, 2. The set T\y spans in this case the whole Hilbert space A(C 2K tg> 
C 2K ), ergo W is optimal. 
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5.2.2 Fermionic Slater witnesses and positive maps 



It is interesting also to consider positive maps associated with fermionic Slater witnesses via 
the Jamiolkowski isomorphism [18]. Such maps employ W acting in 1~La®TLb = C 2K ®C 2K , and 
transform a state p acting in Ha 8> Tic = C 2K (g> C 2K into another state acting in Hb§ Tic = 
C 2K ® C 2K , via the relation M(p) = Tr A (Wp TA ). Obviously, such maps are positive on 
separable states: When p is separable, then for any \$>) G Tis <8> Tic, the mean value of 
(*S>\M(p)\*$>) , becomes a convex sum of mean values of W in some product states \e,f) £ 
TLa <S> Hg. Since VF acts in fact in the antisymmetric space, we can antisymmetrized these 
states, i.e. |e, /) — > (|e, f) — \f, e)). Such antisymmetric states have, however, Slater rank 1, and 
all S1W of class k > 2 have thus positive mean value in those states. This class of positive maps 
is quite different from the ones considered in Refs. [|64| , |65|] ; they provide thus an interesting 
class of necessary separability conditions. The particular map associated with the witness ( p0| ) 
is, however, decomposable, i.e. it is a sum of a completely positive map and another completely 
positive map composed with the transposition. This follows from the fact that W^ A = Q is 
a positive operator and from M(p) = Ti A (Wp TA ) = Tr A (Qp) = M Q (p TA ) = (Mq o T A )(p) 
where M Q (X) = Tt a {QX Ta ) is a completely positive map because Q is a positive operator. 



6 Quantum correlations in mixed bosonic states 

Here we present new results that generalize the previous section to the case of bosons. The case 
of a two-boson system with a two-dimensional single-particle space, where the Slater number 
of a given mixed state can be calculated directly, has already been discussed (c.f. section ||). 
The concern of this section is thus a single particle Hilbert space C K with K > 2. We discuss 
also the relation of the Slater number and positive definiteness of the partial transpose. 



6.1 Slater number of mixed bosonic states 

Following the results from the previous section we may define formally the Slater number k 
of a bosonic mixed state. 

Definition |5.ip . Consider a density matrix p of a two-boson system, and all its possible 
convex decompositions in terms of pure states, i.e. p = YliPiWl 1 (V'lo where ri denotes 
the Slater rank of 1^)1* ; the Slater number of p, k, is defined as k = min{r max }, where 
''max is the maximum Slater rank within a decomposition, and the minimum is taken over all 
decompositions. 

Most of the results concerning Slater numbers for fermions can be transferred directly to 
the Slater number for bosons. If we denote the whole space of density matrices in S{C 2K (&C 2K ) 
by Six, and the set of density matrices that have Slater number k or less, by then Slk 
as in previous section is a convex compact subset of Sljc; a state from Slk will be called a 
state of (Slater) class k. Sets of increasing Slater number are embedded into each other, i.e. 
Sli C SI2 C ...Slk--- C SIk- In particular, Sl\ is the set of states that can be written as a 
convex combination of elementary Slater permanents. There exists, however, one very basic 
difference between bosonic and fermionic systems; in the latter case, namely, the concept 
of partial transposition seems to be completely inappropriate. On the other hand, since 
in bosonic systems Sl\ is the set of separable states, partial transposition provides here a 
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necessary condition for a state to belong to Sl\. For bosonic systems it is thus reasonable to 
ask questions of the sort: Do bosonic entangled states with positive partial transpose (PPTES) 
exist? 

For the case of K = 2 a necessary and sufficient condition for a mixed state to have a 
Slater number one is given either by the Peres-Horodecki criterion, i.e. by the positivity of 
the partial transpose, or by the Wootters approach as discussed in section || The latter result 
we state again in the following theorem: 

Theorem |5.1] '. Let the bosonic mixed state acting in S(C 2 (g) C 2 ) have a spectral decompo- 
sition p = Yli=i where r is the rank of p, and the eigenvectors belonging to 
nonzero eigenvalues Aj are normalized as {9i\^fj) = \8ij. Let = Yl a b v l^ a %\£L) in some 
basis, and define the complex symmetric r x r matrix C by 

C ij = J2^ M <A (83) 

abed 

(e 01 = — e 10 = 1, and e aa = otherwise) which can be represented using a unitary matrix as 
C = UCdU T , with Cd = diag[ci,C2, ... ,c r ] diagonal and \c\\ > | C2 1 > ••• > \c r \. The state p 
has Slater number 1 if and only if 

r 

M<^N- ( 84 ) 

i=2 



6.2 Bosonic Slater witnesses 

We now turn to the more difficult case of K > 2 where as for fermions no complete char- 
acterization of pure and mixed state quantum correlations can be given. Again the concept 
of witnesses turns out to be useful here and we can directly translate most of the fermionic 



results from section ^2. We will only list the theorems here and skip the proofs which are 
discussed in detail in 



Proposition 5.1'. Any bosonic state of class k, p^, can be written as a convex combination 



of a density matrix of class k — 1 and a so-called bosonic A;— edge state 5: 

Pk = (1 -p)p k -i +pS, 1>P>0, (85) 

where the edge state 5 has Slater number > k. 

This again follows from the convexity and compactness of the sets Slk- As in the case of 
fermions, there is an optimal decomposition of the form ((8^) with p minimal. Alternatively, 
restricting ourselves to decompositions pk = YliPi\^Yi with all < k, we can always 
find a decomposition of the form ( pB] ) with 5 £ Slk- This allows to define A;-edge states and 
fc-Slater witnesses as in fermionic case. 



Definition 5^2'. A bosonic fe-edge state 5 is a bosonic state such that 5 — e\ip <k )(tp <k \ is not 



positive, for any e > and \ip <k ). 



Criterion 5.1'. A bosonic mixed state 6 is a bosonic fe-edge state if and only if there exists 



no \il> <k ) such that \ip <k ) G R(S). 
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Definition |5.3p . A hermitian operator W is a bosonic Slater witness of class k (A;-S1W, k > 2) 
if and only if Tr(Wa) > for all a G Slk-i, and there exists at least one bosonic p S S7fc such 
that Tr(Wp) < 0. 



With the same definitions as in the paragraph before lemma 5A in section 52 (but now 
all the operators act in the symmetric space S{C K ®C K )) we have 

Lemma |5.1] '. Given a bosonic /c-edge state 5, then 

W = P--C (86) 
c 

is a bosonic fe-SlW which detects 8. 

Similarly as in section [T^, the simplest choice of P and C consists in taking projections 
onto K(5) and the identity operator on the symmetric space l s , respectively. 



Proposition 5.2'. Any bosonic Slater witness can be written in the canonical form: 

W = W -el s , (87) 
such that RiW) = K(6), where 5 is a bosonic k-edge state and < e < inf|^<fc^ (^'|W'|^') ■ 

6.2.1 Optimal bosonic Slater witnesses 

We will now discuss the optimization of bosonic Slater witnesses, applying the same definitions 
of finer and optimal as for fermions. Again we have that if the set Tw of pure tangent states 
of Slater rank < k, 

T w = {|V?) such that (V>|W|V>> = 0} , (88) 

spans the complete Hilbert space S(C K (8>C ), then W is an optimal A>S1W. Otherwise we can 
optimize the witness by subtracting from it a positive operator P, such that PT\y = 0. For ex- 

wn 1/2 

ample, for Slater witnesses of class 2 this is possible provided that inf| e ^ gC K [P e W e P e ] m in > 
0. Here for any X acting on S{C K C K ) we define 

(e,.\X\e,.) + (e,.\X\.,e) 
+ (.,e|X|e,.) + (.,e|X|.,e) 

as an operator acting in C K , and [AT] m i n denotes its minimal eigenvalue (compare |34)). An 
example of an optimal bosonic witness of Slater number k in S{C K (g> C K ) is given by 

W = l s - -j^P , (90) 

where V is a projector onto a "maximally correlated state", \^>) = -7= X^Li &o&o|^)- That 
this witness is optimal can be seen as in section p.2.1 . 



Finally we only want to mention that it is again possible to relate positive maps to bosonic 
Slater witnesses, leading to another interesting class of necessary separability conditions. As 
for fermions (c.f. section 5.2. 1| ) the map associated with the witness (90) is decomposable. 
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6.3 Bosonic entangled states and positive partial transpose 

In the case of distinguishable particles there are strong indications that the Schmidt rank 
of entangled states with the PPT property is not maximal, and rather low. This has been 



conjectured for the states in a 3 x 3-dimensional space in Ref. [60|, and for the states in a 



N x A-dimensional space in Ref. [67]. We expect that similar situations take place for the 
Slater rank of bosonic states. In fact we have not been able to find a two-boson mixed state 
that would have the PPT property and would be nevertheless entangled, i.e. it would have 
Slater rank strictly greater than 1. For bosonic 3-qubit states as well as for bosonic iV-qubit 
states of rank < N we were even able to prove that such states do not exist. 

Our findings are summarized in the following two theorems. The first one deals with the 
case of K = 3, i.e. three-dimensional single-particle space, but can be easily generalized to 
arbitrary K. The second one deals with bosonic systems of N qubits, i.e. N bosons in a 
two-dimensional single particle space. 

Theorem 6.1. Let a bosonic mixed state p acting in S(C 3 & C 3 ), which we treat as a sym- 
metrized space of Alice and Bob, have a positive partial transpose with respect to, say, Alice. 
Then 

I. if the rank of p, r(p) = 3, then the Slater rank of p, Sl(p) = 1, i.e. p is separable. 

II. if the rank of p, r{p) = 4, then generically there exist exactly 4 product vectors of the 
form \e{,ei) in the range of p, and the Slater rank of p, Sl(p) = 1, i.e. p is separable, 
provided (ej, ej|p _1 |ej, ej) = Pi5ij, where the inverse is taken on the range. In such case 

P — /Zj=l Pi l^i) ^i) (^*> I ■ 



Proof. This theorem is a simple generalization of the result of Ref. | ]68| . In particular from 
Theorem 1 of that paper follows that any (symmetric or not) PPT state of rank 3 is separable, 
which proves |]. To prove we observe that the vector |e,e) is in the range of p if it is 
orthogonal to vectors which span the kernel of p. There are two such vectors \cj>j) , j = 1,2, since 
the whole space 5(C 3 (8> C 3 ) is 6 dimensional. Since in a convenient basis up to normalization 
we can write |e) = (l,zi,z%), the equations (4>j\e, e) = constitute a system of two quadratic 
polynomials in z\ and z-i- Such a system of equations has generically 4 solutions |ej, ej), 
i = 1, . . . ,4. From that it follows that p is separable if and only if p = Yli=iPi\ e ii e i){ e i> e %\- 
Using this equation, we observe that if p~ 1 \ei, ej) = \v), then |ej, ej) = p\v) and from the linear 
independence of |ej,ej) we get that (ej,ej\v) = for j ^ i, and (ei,ei\v) = 1. This implies 
that (e^ ei|p _1 |ej, e^) = Pi8%j. 

Theorem 6.2. Let a bosonic mixed state p acting in the symmetric A-qubit space S(C 2 (8> 
. . . (8> C 2 ), which we treat as a symmetrized space of Alice, Bob, Charlie etc., have a positive 
partial transpose with respect to, say, Alice. For N > 3 let r(p) < N. Then p is separable. 

Proof. Let us consider first the case N = 3. The state p is bosonic, i.e. has a maximal 
rank 4, which is the dimension of the space of totally symmetric vectors. Suppose that the 
state p is, however, also a PPT state with respect to Alice. We can treat it then as a PPT 
state in a 2 x 4 dimensional space of Alice and Bob and Charlie grouped together. Such a 
state of rank 4 is necessarily separable and can be represented as p = ^2 i= iPi\ei,^>i)(ei,^i\, 
where the projectors onto the vectors are in general entangled pure states of Bob and 
Charlie Since they belong to the range of the reduced density matrix Tr^/3, they must 
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belong to S{C 2 <S> C 2 ). Obviously, the vectors |ej,\l/j) belong to the range of p, and thus to 
the symmetric space S(C 2 (g) C 2 (8> C 2 ). Direct inspection shows that the only vectors of this 
form are |ej,ej,ej): writing the Schmidt decomposition = a|00) + b\ll), and expanding 
|e,) = q|0) + /3|1) in the same basis, we immediately see that both ab and /3a must be equal 
to zero. 

The proof for higher N is analogous. We assume that p has rank < N and a positive 
partial transpose with respect to Alice. Since p, regarded as a matrix acting in a 2 x 2 N ~ l - 
dimensional space, is supported in reality in C 2 ® ^((C 2 )®^ -1 ), i.e. in a 2 x A-dimensional 
space, the results of Ref. |6^] hold and p must be separable with respect to the partition Alice 
- the rest, and therefore have a form p = Ylf=iPi\ e it ^i)( e «; But then, using the analogous 
argumentation as above we see that | , Vl/j) = e», . . . , ei). 



7 Conclusions 

In conclusion, we have studied and characterized quantum correlations in systems of (i) two 
distinguishable particles, (ii) two fermions and (hi) two bosons. As an equivalent to the 
well-known Schmidt decomposition and rank for distinguishable particles we introduced the 
concepts of Slater decomposition and rank, i.e. writing fermionic and bosonic pure states as 
sums of orthogonal Slater determinants or permanents, respectively. 

For the case of the smallest non-trivial Hilbert space we were able to identify interesting 
properties shared by all these systems: an appropriate time reversal operator allows to define 
a correlation measure, the concurrence, which can also be expressed through the coefficients 
of a state in the so-called magic basis. This also allows for perfect analogues of Wootters' 
formula \27\ for the correlation measure of mixed two-qubit states. As we noted, these analogies 



follow directly if the system of two fermions in a four-dimensional single-particle space is 
considered as the most general one from which the other two can be derived. 

Subsequently we analyzed systems of fermions or bosons in higher-dimensional single- 
particle spaces and systems of more than two indistinguishable particles. There is no straight- 
forward extension of the results obtained for low-dimensional Hilbert spaces. We gave nec- 
essary and sufficient criteria to identify the Slater rank of pure two-particles states by con- 
tractions of the e-tensor with its coefficient matrix. These results are sufficient to classify 
quantum correlations of pure states in arbitrary-dimensional single-particle spaces. For the 
case of three or more indistinguishable particles it was first established that a general state 
cannot be written as a sum of orthogonal Slater determinants, which demonstrates yet an- 
other analogy to distinguishable particles. We showed that non-correlated pure states can be 
identified through projections of the multi-particle state onto states of two fermions. 

For general two-fermion mixed states we reviewed the concept of witnesses of Slater number 
k and introduced a similar concept for bosons. In both cases we showed how to optimize these 
witnesses and gave an example of an optimal witness. 

Finally we explored the role of the PPT criterion in bosonic systems. We gave separability 
criteria for low-rank two-boson mixed states in a three-dimensional single-particle space having 
a positive partial transpose. We also showed that there are no bosonic three-qubit states that 
have a positive partial transpose and nevertheless contain quantum correlations. The same is 
true for bosonic A-qubit states of rank < N. 
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A Mode correlations 

Here we discuss briefly the rather complementary concept to deal with fermionic correlations 
introduced by Zanardi [^l| . Again consider the quantum dot example and suppose the system 
evolves to a state (here we neglect interaction of the particles) 

W) = J[|^>IT>®l^>U> + lx)IT>®lx)U> + l^>U>®l^>|T> + lx>U>®lx>IT>] (91) 



2 
1 



|0) (92) 



V2 

If the tunneling barrier is raised in this state then either Alice or Bob find to electrons in their 
dot, i.e. there is no particle entanglement present here. But if we forget about the spin and 
map the two particle basis state ///j|0) onto \n^) (g) \n x ), where and n x are the occupation 
numbers of the left and the right dot, respectively, then the result is a state oc |2)<g>|0) + |0)<8>|2). 
This is an entangled state, but now the entangled entities are not the particles themselves but 
"modes". Each mode is associated with (in this case more than one) creation operator fj. 

This is the concept of mode entanglement. In pl| fi"2l correlations are analyzed not by 
looking at the particles themselves but by partitioning the state space into different modes. 
In the simplest case these modes are associated with the 2K single particle basis states such 
that the occupation number is either or 1 in the case of fermions. These 2K modes are then 
mapped onto the space of 2K qubits such that the k th qubit is in state |1) (state |0)) if the 
k th mode is occupied (not occupied). 

More formally this mapping is describe for the basis states of N fermions as [41] 

A: H£ K (N) -► (C 2 )® 2K (93) 

ft---fl N \0) - ®llx\x(k;h,...,i N )) (94) 

(for i\ < ... < in) and linearly extended to the complete Hilbert space. The characteristic 
function x ls defined as 

x(*;«i>-.«*) = ( o k^{h,...,i N } (95) 

The mapping is not bijective in this form as not all possible i'T-qubit states are realized, 
i.e. ?i2 K (N) is not isomorphic to this space. An isomorphism can be formed by not fixing the 
number of fermions in the system, i.e. by considering 

A;H£ K ^ (C 2 f 2K (96) 

where H$ K = 0^L 1 W^(A r ) is the total Fock space of 0, . . . , 2K fermions. Notice that this 
point might raise problems in concrete applications as it requires sources and drains to adjust 
the number of particles. 



35 



B Proofs (Quantum correlations of pure fermionic and bosonic 
states in higher-dimensional Hilbert spaces) 

Here we give proofs of some of the lemma on indistinguishable particles in higher-dimensional 
Hilbert spaces from section 

B.l Fermionic states 

First notice that according to equation (^) for any complex antisymmetric 2K x 2K coefficient 
matrix to there is a unitary transformation IA such that 

^. = (UwU T ) {j = Zl ^ M ~W + 22E 12ijM...(2K) + + Zk£ 12..,j (9?) 

where e ll >— * 2Jf is the totally antisymmetric unit tensor in 2K x . . . x 2K and z k > 0. 

Proof of lemma |4.1| : [=>] Consider the case a = 5, (3 = 6. In the sum all terms with at 
least one i,j,k,l £ {5,6} vanish, and thus 

6 4 

mjw kl e^ km = ^kie llkl = C (|t;)) (98) 
i,j,k,l ij,k,l 

where |t,) = VijfJfjlO) e A{Hi (8) 7U) and Vij is constructed from t«y by deleting the last two 
columns and rows. 

Because by assumption \w) has Slater rank 1 also \v) has to have Slater rank 1 unless 
\v) = 0. In both cases C (\v)) has to vanish. The same argument can be repeated for every 
combination of a and (3 (only deleting different rows and columns from the coefficient matrix 
w). 

[<=] We will show the negated claim: there is at least one sum that does not vanish if \w) has 
Slater rank > 1. The two possible cases to be distinguished are Slater rank three or Slater 
rank two. If the Slater rank is three then w has full rank and thus det(w) ^ 0. Using the 
expression for the determinant of an antisymmetric matrix from equation (|l8|) it holds 

6 6 

48detH= W <*P E u>ijWkie ijklaf} + 0, (99) 

a,/3=l i,j,k,l=l 

y v ' 

x a/3 

implying that at least one x a/3 ^ 0. 

If | to) has Slater rank 2 then r{w} = 4 and we can find a unitary transformation U such 
that 

6 

w'i'ji = ^2 U i'i w ij U j'j ( 10 °) 

contains only zeros in the last two rows and columns. \w) has Slater rank 2 if and only if 
\w') = Y^ij w ij(f')\ (/')!• |0) nas Slater rank 2 and it follows (to' is a 4 x 4 matrix only and thus 
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it is possible to use the concurrence from section |3] which does not vanish if and only if \w') 
has Slater rank 2) 

4 

O^C(K)) = E ^'mn^™ (101) 
i ,j,m,n=l 

= E <^n^ mn56 (102) 
i ,j,m,n=l 
6 

= E Uui WiijiUjjiU mm iw m i n iU nn i e^ mn5G (103) 

i j',m,?i=l 
i' ,j' ,m' ,n'=l 

6 

= E W raV ^ii' ^j'j' ^mm' ^nn' £ 4Jimn56 (104) 

i,j,m,n=l _ r~ 

p i j m n 

i' ,j' ,m' ,n'=l 

It is easy to verify that £ l i mn \ s a l so antisymmetric in all indices so it has to have the form 

p'j'm'n' = £ 56 /j'mV56 + ^/j'm'n'46 + _ _ _ + ^i'j'm'n' 12 

(where some but not all £ a/ 3 may be zero) and then 

6 

n -L \ „ , „„ (t Ji'i'm'n'56 , /- A'j'm'n'46 . , ,j A'j'm'n'12 \ 

7^ 2^ w i'j' w m'n' ( ?56£ +^46^ + ■ • • + 4l2£ ) (106) 

i' ,j' ,m' ,n'=l 
6 6 

^ E^ E ^ lt ^w vf w m , (107) 

a/3=l i',j',k',l'=l 

S v ' 

such that again at least one x a P ^ 0. 



Proof of lemma ^2 : li N = K then the claim is that |u>) does not have maximal Slater 
rank if and only if 

2K 

E ^•••^ 2 k-i^ 1 - 42K =0 (108) 
h,— =1 



This is the statement of equation (18) so lemma 4.2 surely holds for all cases [K = N,N]. This 



is the starting point for giving the proof by induction. For a combination [K, N] (N < K) 
the assumption will be that the lemma holds for [K, N + 1] , [K — 1, N] and [N, N] (figure |H| 
illustrates the induction process). 

Introduce a two-fermion state \v) £ A{'H.2K-2® r H-2K~2) by constructing its coefficient 
matrix v from w by deleting the rows and columns a 2 (K-N) an< i a 2(K-N)-i- Since is as- 
sumed to have Slater rank < iV also \v) has to have Slater rank < N (because its mathematical 



37 



K 

2 3 4 5 6 



2 
3 
4 



N 



Figure 5: Construction of the proof of lemma 4.2 by induction. Black dots mark the combi- 



nations [K=N,N] where equation (^) can be applied. 



rank cannot raise when some rows and columns are deleted). Then it follows 

2K 

]T w^.-.w^^e^-^^-^-^ = (109) 

h,— :*2JV=1 

2(K-1) 

^^■■^- lis / ,1 '- !2 " ai - a21 ' f - 1 - w) = o, (110) 

ill— i*2iV = l 

because by assumption the lemma holds for K — 1 and N. 

[<=] We will show that there is at least one non-vanishing sum if the Slater rank is > N. 
Then if all terms vanish the Slater rank has to be < N. Again there are two cases to be 
distinguished: 

Case 1: \w) has Slater rank > N. By assumption the lemma holds for [K,N + 1] and 
thus 

2K 

+ £ w llt2 . . . ™ i2N+lW ^--^ +2 ^--- 2 (K-;v-i) (m) 

il,„. ,J2JV+2 = 1 

for at least one combination of a%, ... ,ct2(K—N-l)- But then 

2K 2K 

fl -/ 111- ID- ■ 111- p U---«2JV*2JV+l«2JV+2ai---02(A'-JV-l) Hi 91 

U T 2-^1 W l2N+H2N+2 2-~t 1112 ' ' ' W *2iV-l«2JV fc \ 11Z ') 

«2JV+l,i2AT+2 'b-.'2JV 

N v ' 

x i 2N+l i 2N+2 a l--- a 2(K-N-l) 

such that at least one x i2N+li2N+2ai - a2 ( K - N ~^ / 0. 

Case 2: \w) has Slater rank N. Then r{w} = 2N and there is a unitary transformation 
U € SU(2K) such that 

w'irf = £ Ui'iWijUj'j (113) 
ij=l 
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contains zeros in the last (2K — 2N) columns and rows. \w'} also has to have Slater rank N 
and because the lemma is assumed to be valid for [K = N, N] it holds 



2N 



^ w'- ■ w' ■ e h - i2N (114) 

i*2JV = l 

2K 

V w' ■ w'- £ ii-i 2N (2N+l)...(2K) (U5) 

ill— ,*2JV = 1 
2K 

E ^«i 2 ^--^^^ (116) 

i\...i 2N = l 

ii---i 2 jv =1 

2K 2K 

= y w fr ...w r y u hjl ...u i2Nf £ ii-iM™+i)...(2K) (117) 

i. ■< J\J2 J2N-1J2N ' ■> l 2Nj 2N v ' 

ji-j' 2N =l h...i 2N =l 

-][... j' 2N {2N+l)...(2K) 

gj'i---j' 2N (' 2N + 1 )---( 2K ) h a s ^ antisymmetric in all indices so it has the form 

2K 

~f 1 ... j ' 2N (2N + l)...(2K) = £ ^...^(K-iV)^-^" 1 - 02 ^^' (H8) 

"ll— , a 2(K-N) =1 

(some but not all £ might be zero) and then 

2K 2K 

M y ^... a2{K _ N) y "^•••«:^ v. (119) 

ai-a 2 (x-iV)=l j;...i 2JV =l 

V V ' 

a . a l--- a 2 (A'-iV) 

where at least one x ai "' a2< - K - N ^ / is necessary. 

Proof of lemma |4.3| : [=^] \w) has Slater rank one so there exists U, UU^ = 1, such that 

w' ijk = w lmn U u U mj U nk = Z£ vMm...2K (12Q) 

Then, for arbitrary |a) and writing \a'} = U^a), we have 

2K 2K 

\v') = R W) \w') with ^ = 3 y w' ijk a' k = y ze^-^aVh = {A® 0)y (121) 

k=l k,l=l 

with the 3x3 antisymmetric matrix A and as the (2K — 3) x (2K — 3)-dimensional zero. A 
has rank zero or two and thus there exists unitary V such that 

{Vv'V T ) l3 = /?^' 345 -( 2 *) (122) 
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for some (3. Inserting equations ( |12C| ) and (121) shows that \v) = Ri a \\w) indeed has Slater 
rank one (if (3 ^ 0): 

2K 

0£ ij34 5 ...(2K) = V im w opq U mo U np U kqai U^V jn (123) 

l,m,n,o,p,q=l 
2K 

= £ ( uTy )oiWo Pq a q {U T V) pj . (124) 

o,p,q=l 

[<£=] Here we make use of the fact that an iV-fermion state \w) (coefficient matrix w) in a 2K- 
dimensional single-particle space has Slater rank one (it can, after some transformation of the 
single-particle basis, be written as a single Slater determinant) if and only if there exists a set 
of orthonormal vectors {e a } OL= i K with eP-e 8 = 5 a/3 , such that for a\, . . . , on G {1, . . . , 2K} 

w ijk ef 1 ... e a k N = z e "i-«;v(iv+l)...(2*0 ( 125 ) 

These vectors are the columns of the matrix U that transforms the single particle space. 

Choose an arbitrary e 1 , \e l \ = 1, such that Vj k = Yli w ijk e i is not the 2-fT x 2-fT-dimensional 
zero. The resulting two fermion state \v) by assumption has Slater rank 1 and thus there exist 
e 2 , e 3 , . . . , e 2K such that e. a e^ = 5 a/3 and 

2K 2K 

£ = £ = (126) 

j,k=l i,j,k=l 

But we furthermore claim that \w) has Slater rank one: 

2K 

£ WijkefefeZ = ze a ^ 4 -^ for a, (3, 7 G {1 . . . 2K} . (127) 

i,j,k=l 



Because of equation ( |126j ) this is obvious if (a, /3, 7) is a permutation of (1, 2, 3) and otherwise 
if a = 1. So further attention has only to be paid to the cases where 1 < a < (3 < 7. 

Consider with = Yl w ijk e k for 1 < a < 3. By assumption the two-fermion state 
\v a ) has Slater rank one (or is zero) and thus there is at most a two-dimensional subspace A4 
such that for any two a G M. and b G M., a 7^ b, Yl w ijk^t a j^k 7^ 0- Because of eqn. ( |126j ) 
M. is spanned by , e 7 where {/?, 7} = {1,2,3} \ {a}. The e's are orthogonal and thus 
Wijke^ej e\ = 0ifl<a<4 and «</?'< 7'. 

For 3 < a < (3 < 7 fixed suppose ^ Wy&ef 7^ 0. Then there exist a two-dimensional 
subspace spanned by e 2 and e 3 and an orthogonal subspace M spanned by e@ and e 7 
such that Yl w ijk [e/ + e"] ajb k / for a and b both in A4 or both in M and a ^ b. This 
contradicts the assumption that i?| e i^ + | eQ ^ has Slater rank one. 

Proof of lemma |4.5| : [=^] \w) has Slater rank one, i.e. we can find e 1 , . . . , e 2K such that 

2K 

£ v Hl ... iN e£...eff = ze*-*«l N+ Q-<™ (128) 
ill— >»jv=1 
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Now consider fj 1 ...j iV _ 1 = ^ u>j 1 ...j Ar Xj iV for arbitrary x. If x = e a then \v) has Slater rank 1 if 
a < N or is zero if a > N. Otherwise it holds that 

2K 2K 

E %...i JV _ie i "_ 1 = fora>iV and E ^iiv-i = ( 129 ) 

*jv-i=1 ijv-i=l 

and we can choose 6 a = e a for a > N, b N = x and the remaining 6 a for a < N such that 
6«i6«a = fiction tQ get 

2a: 

E ^..^- 1 C-"C" =^^- lN " W (130) 

ii,... ,ijv-i=1 

which means that |u) has Slater rank one if z ^ and is zero otherwise. 
[<=]: similar to the three- fermion case. 



B.2 Bosonic states 



Proof of lemma 4.1: [=£-] Let \v) have Slater rank one and for any 1 < a < 3 let v' be the 
2 x 2-matrix that results from deleting the row and column a. If \v) has Slater rank one then 
also |u') = Y2ij=i v ijKb\\0) has Slater rank one (unless it is zero) and the claim follows from 
equation (|30"|). 

[<=] We will show that at least one sum in equation (^) does not vanish if \v) has Slater rank 
larger than one. The two possible cases are Slater rank three and Slater rank two. If the 
Slater Rank is three then v has full rank and det v ^ 0. There exists a unitary transformation 
U of the single particle space such that 

3 

E Uu'Vi'j'Ujji (131) 
i>i=i 



V 'i'j' 



is diagonal. In this basis 

3 3 

0^3!deW = E4 E ^^"^ ( 132 ) 

i=l fc,m=l 
3 3 

= E^ E U kk >v k >i>U k i>U mm >v m i n 'U mn 'e lkm e tkm (133) 

i=l fc,m=l 

fc',i',m',n'=l 

3 3 3 

= E^E^4 E v ™ Vm'n'S ak ' m 'e 0n ' ^ (134) 

i=l a,/3=l k' ,1' ,m' ,n' =1 



and thus at least one ^ 0. To show that his implies that at least one x aa ^ 0, let 
{a, /3} S {1,2,3} and let 7 be the remaining index. Then 

| x a/9| = 2|e Q ^^^ Q v 77 + e a ^e /37a w /37 i; 7Q | (135) 

= 2|e^e /3 ^v /3Q v 77 -e a/37 e /3a7 w / 3 7 v 7Q | (136) 

1 v 2 

= 21^^-^^0,1=0 44> 1 = ^77 2 ^ ' ( 137 ) 
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On the other hand if a £ 1, 2, 3 and /?, 7 are the remaining two indices, then 



(138) 
(139) 



and thus if x 11 = 0, x 22 = and x 33 = then v 2 s = ^22^33, ^13 = ^11^33 and v\ 2 = v \\V22 ■ 
Then equation (137) is fulfilled and x a ^ = for all choices of a and 0. Thus at least one 
x a/3 q i m pij es that at least one x aa 7^ 0. 

If the Slater Rank is two then there exists a unitary transformation U such that 



^2 Uu'Vi'jrUjf = diagCv'^v'^O), 



and thus 



Jikm 



^ v 'kk v 'mm £3km £ ' 

3 

ET T T T ~T~ T T T 3 AjTTI / ' lit 

<Jkk' v k'l' u kl> u mm' v m'n' u mn> e e 



k,m=l 
3 



(140) 

(141) 
(142) 



k,m=l 



such that the above arguments can be repeated. Using this, lemma f4.2| can then be proven by 
induction. 



Proof of lemma 4^3: We will give the proof for the case N = 3. [=£-]: \v) has Slater rank 



one and, following definition 4.1, there exists a unitary transformation U such that 

K 

J'ijk 



v'nh = ^ v lmnUilUj m U kn = zSixSjiSkl- (143) 



Lm,n=l 



Then for arbitrary \a) it is obvious that i?i a \|u) has Slater rank one (unless it is zero), because 
(here a' k = U M ai) 



K 



3 Yl V 'i3k a 'k = a l z ( 144 ) 



k=l 



[4=] Notice that if a K x K-matrix v has rank one and (normalized) eigenvector e, then if 
the normalized vector x maximizes || vx || then x = e. Now let x and y be normalized 
vectors varying over the unit sphere and maximize || Ylijk v ijk x iyj II- A maximum exists 
because the product of two unit spheres is a compact space and || . || is a continuous function. 
Let the maximum occur at x = and y. Now y is the eigenvector of Yli v ijk e i > which 
by assumption has rank one, because y maximizes || Yli j v ijk^yj II- Similarly e^ 1 ) is the 
eigenvector of ]T\ Vij^yi- Then we have 



^Vijkef^yj = XxVk and y^^ijkUief 3 = Aiej^ (145) 
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and by comparing both equations 

3 



y = e^, Ai = Ai, ]T v^efUfe^ = A 1; (146) 

where Ai can be chosen real and positive by adjusting the phase of e*- 1 ). Furthermore 

Yli,j v ijk4 Zj = for any z 1 e^ 1 ). 

Now repeat this maximization, but this time restrict x and y to the space perpendicular 
to e^ 1 ). If the result is called e^ 2 \ then Yli j k v ijk e i^ e ^ e k^ = A2 and A2 < Ai. In this way a 
complete basis is constructed. If v is written in this basis then 

= feW + ... + eW) = -Ldiag[Ai,... , X K ] , Ai > A 2 > . . . > X K > (147) 



By assumption has to have Slater rank one, i.e. v' has to have rank one, and it follows 
that A2 = ■ • • = Xk = 0. Thus 

E v m4 a) ef } ei 7) = XiSxaSysS^ (148) 
and |u) has indeed Slater Rank 1. 
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